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ABSTRACT 



The surface flux transport (SFT) model of solar magnetic fields involves empirically well-constrained velocity and magnetic fields. 
The basic evolution of the Sun's large-scale surface magnetic field is well described by this model. The azimuthally averaged evolution 
of the SFT model can be compared to the surface evolution of the flux transport dynamo (FTD), and the evolution of the SFT model 
can be used to constrain several near-surface properties of the FTD model. 

We compared the results of the FTD model with different upper boundary conditions and diffusivity profiles against the results of the 
SFT model. Among the ingredients of the FTD model, downward pumping of magnetic flux, related to a positive diffusivity gradient, 
has a significant effect in slowing down the diffusive radial transport of magnetic flux through the solar surface. Provided the pumping 
was strong enough to give rise to a downflow of a magnetic Reynolds number of 5 in the near-surface boundary layer, the FTD using 
a vertical boundary condition matches the SFT model based on the average velocities above the boundary layer. The FTD model with 
a potential field were unable to match the SFT results. 
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1. Introduction 

The flux transport dynamo model (FTD) attempts to explain 
the large-scale evolution of the Sun's magnetic field. The cen- 
tral ideas behind the model are that poloidal flux is wound up 
by differential rotation until it becomes sufficiently strong that 
magnetic buoyant flux tubes emerge through the solar surface. 
The erupted field is in the form of a bipolar active region, and 
the two opposite polarities are observed to be systematically 
tilted with respect to the equator (Joy's law). This tilt is such 
that the leading polarity is slightly closer to the equator than 
the following polarity. This latitudinal offset means that poloidal 
field has been created from the toroidal flux and, in the lan- 
guage of dy namo theory, the emergence pr ocess is a non-local 
alpha effect (iKitchatinov & O lemskov 201 lal) . The poloidal flux 
is then stretched and diffused by surface motions, reversing the 
polar fields and completing one half of a solar cy cle. For a re- 
view of the basic ideas, see ICharbonneaul (|2010|) . This picture 
has recently gained obse rvational sup port from the analysis by 
Dasi-Esp uig et al.l (T2010 ) and later by IKitchatinov & Olemskovl 
(2011b), which show that the observed sunspot group tilt an- 
gles, which go into the construction of the poloidal source term, 
vary systematically from cycle to cycle in a way which possibly 
can explain the observed changes in cycle amplitudes during the 
twentieth century. 

The winding up of the field by differential rotation and the 
rise of the tubes to the surface are hidden below the photosphere. 
The evolution of the field after it has broken through the sur- 
face can be and has been observed. The surface flux transport 
model (SFT) has been found to provide a good description of the 
large-scale e volution after e mergence. For a detailed historical 
account, see ISheelevI (120051) . This model assumes that the mag- 
netic field is purely radial at the surface and evolves passively 



driven by surface flows including differential rotation, merid- 
ional circulation and small-scale convective motions (granula- 
tion and supergranulation). The small-scale motions essentially 
cause the magnetic field to undergo a random walk and hence 
can be treated as a diffusive term. The ingredients which go into 
the SFT model are all observable, as is the output of the model 
- it is thus tig htly constrained an d supported by observations. 
For example. ICameron et al.l (12010b showed that the SFT model, 
with the observed cycle-to-cycle variations of the tilt angle, can 
reproduce the inferred open magnetic flux of the Sun. Since the 
open flux during the maxima and minima of activity reflect the 
equatorial and axial dipole moments respectively, the model's 
ability to reproduce the open flux over an extended period is a 
strong test of the model. 

In this paper we investigate what constraints can be inferred 
for the FTD model, given that it should also reproduce the same 
surface dynamics as is described by the SFT model. We have 
used the FTD code developed at the Max-Planck-Institut fiir 
Sonnensystemforschung. For the SFT model we have used a 1-D 
surface flux transport model developed at the MPS. The 1-D SFT 
model includes exactly the component which can be compared 
between the two models. The details of the two approaches will 
be discussed in Sect. 2. In Sect. 3 we present the results of the 
simulations and compare the two models. The effect of varying 
some of the most important unconstrained parameters and the 
boundary condition will be discussed in Sect. 4. We will con- 
clude in Sect. 5 with the finding that the appropriate boundary 
condition for FTD models is that the field is vertical at the sur- 
face, and that a certain amount of turbulent pumping must be 
included for the FTD simulations to mimic the surface behavior 
of the SFT model and to thus match the observations. 

Downward pumping has in particular been discussed for 
the base of the solar convection zone. Direct numerical simu- 
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lations show a downward transport of large-scale magnetic field 
near the base of convective unstable layers (e.g..LFennings et al.l 
U992tlTobias et al.lll998ll2001l:IOssendriiver et al.ll2002l) though 
it is not clear whet her this should be interpreted in terms of 
turbulent pumping (Zel dovichlll957t lRadlerlll968l) or of topo- 
logical pumping dDrobvs hevski & Yuf erev|[l974l) . In mean-field 
dynamo models of the solar cycle turbulent pumping is often 
not included. In those cases where it is included, most of the 
attention is focussed on its role in transporting flux into the 
top of the stable layer immediately below the convection zone 
(see for example iBrandenburg et alj 1 1992h iKapylaetaL {2006 ; 
iDo Cao & Brunll201 U iKitchatinov & Olemskovll201 lal |2012|) . 
The effect of pumping throughout the convection zo ne was con- 
sidered by dGuerrero & de Gouveia Dal Pinoll2008l) . who found 
that the pumping affects whether the preferred mode of the solu- 
tion is dipolar or quadrupolar, and identified the possible impor- 
tance of radial transport by the pumping in the dynamo process. 
In this paper we are especially payi ng attention to the pumping 
in the near-surface boundary layer (Miesch & Hindman 2011) 
which is necessary to match FTD to SFT simulations of the mag- 
netic flux on the solar surface. 

2. Physical models and numerical codes 

2.1. The Flux Transport Dynamo (FTD) model 

The flux transport dynamo equations describe the induction, ad- 
vection and diffusion of a large-scale magnetic field. Their ax- 
isymmetric form is: 
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where A(r, 6) is the 0-component of the vector potential associ- 
ated with the poloidal components of B, B(r, 6) is the toroidal 
component of the field, u m (r,ff) is the velocity in the merid- 
ional plane, D(r, ff) is the angular velocity, u p (r, 6) is a veloc- 
ity field corresponding to the pumping of the magnetic field and 
a is a source term in the equation for A corresponding to the 
generation of poloidal flux from toroidal flux. Since the pur- 
pose of the current study is to compare the response of the SFT 
and FTD models to equivalent sources of poloidal flux, we re- 
strict ourselves to the case a = - for other choices of a 
we would need to modify the source term in the SFT model 
accordingl y. In relation to the term « p it is important to note 
that, as in iGuerrero & de Gouveia Dal Pinol (120081) . it does not 
correspond to a true motion of the fluid and need not satisfy 
V ■ pu p = 0. Rather it is a parametrization of the effect of the tur- 
bulent motions on the field: for diamagnetic pumping it has the 
form u p (r,6) = -j^r/. Other effects, such as topological pump- 
ing, are also expected to transport the field downwards, and for 
this study we assume that the combined effects of the turbulent 
convection, including diamagnetic pumping, can be written in 
the form « p (r, 9) = -5V77, with k > 1. This choice allows us to 
vary the magnitude of the pumping in the near surface layers. 



We solve the dynamo equations (Q} and forward in time in 
a spherical shell ro < r < Rq with inner boundary ro = 0.65Rq 
matching to a perfect conductor and outer boundary matching to 
either a radial field or vacuum conditions outside. This leads to 
the boundary conditions 
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for the field to be vertical at the Sun's surface, or alternatively 
A = ^ fl ^(cos0), (5) 
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for matching to a potential field outside. At the poles we require 
regularity resulting in 



B = at 6» = 0,7T 



(8) 



The equations are discretized using second order accurate cen- 
tered finite differences on an equidistant grid and forwarded in 
time with an Alternating Direction Implicit scheme for the dif- 
fusion terms and an explicit scheme for the induction and advec- 
tion terms. The co de is tested against the dynamo benchmark of 
iJouve et al .1 (120081) . 

For current purposes we will consider a = so that there 
is no source of poloidal field during the simulation. From any 
initial condition the field must then eventually decay towards 
zero, however at any finite time the magnetic field will depend 
on the initial field and can be compared with the result of the 
SFT model. 

For the initial condition we take 



5 = 0, 



(9) 
(10) 



where r x = 0.80%. #1 = 80°, 8 2 = 86°, A0 = 2.9° and 
Ar = Q.QIRq. This corresponds to an isolated bipole emerg- 
ing on the solar surface slightly north of the equator. Since both 
the SFT and FTD studied here are linear, the evolution of such 
a bipole is independent of the emergence and evolution of other 
emerging groups. To check that the models are consistent, it is 
thus sufficient to follow the evolution of a single feature starting 
near the equator to the poles. 

The velocity i n the meridional plane is taken from 
Dikp ati et al.1 d2004t) . The velocity components can be written 
in terms of a stream function as 
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For the reference case we take g = 1.5, to = 1.5, p — 3, ro = 0.7 
and frj = £(rn). In all cases uo is chosen so that the maximum 
meridional velocity at r — Rq is 15 m/s. The resulting veloc- 
ity approximates the meridional circulation in the so l ar co nvec- 
tion zone derived by numerica l modeling of Rempel (2005J) and 
iKitchatinov & Riidigerl (120051) and is consistent with the veloc- 
ity in the subsurf ace layers as deriv ed from helioseismology (as 
measured, e.g. bv lGiles et al.l fl997) 

Th e differential rotation is 
(2000), and is also used e.g 



taken from Belvedere et al 
by IKitchatinov & O lemskov 
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where the coefficients cy are given in Table 1 of iBelvedere et al.l 
(2000). This approximates the internal rotation of the Sun as de- 
rived from helios eismological inversions (as reported, e.g., by 
ISchou et alJl998h . 

For the diffusivity we assumed 
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with 770 = 0.1 km s , 771 = 10 km s and 772 = 250 km s , 
see e.g. M unoz-Jaramillo et al ] (1201 11) . Here 772 represents the 
turbulent diffusivity in the near-surface layers, 771 in the bulk of 
the convection zone, and 770 in the overshoot region at the base of 
the convecti on zone. Other choices will be considered in Sect. 4. 

Recently IKitchatinov & Olemskovl d201 lal) have highlighted 
the importance of downward pumping of magnetic fields due to 
gradients in the turbulent diffusivity, and have argued that this is 
particularly important near the base of the convection zone. We 
here consider downward pumping in the near-surface layers. We 
have found it necessary to increase the strength of the downward 
pumping from its usual value of (l/2)(8T]/8r) in order to obtain 
a match between the FTD and SFT models. We have therefore 
introduced in Eqs. (Q} and (0 a scaling factor k which we have 
varied between and 20. The diffusivity profile and the corre- 
sponding diamagnetic pumping velocity with k — 1 is shown in 
Fig. Q] 

While we solve both Eqs. (Q} and (|2]i we note that the com- 
parison with the SFT model only depends on Eq. (Q} as a — 0. 
The physical ingredients which affect A are the meridional flow, 
the radial and latitudinal diffusion, and the downward pumping. 
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Fig. 1. The assumed profile of the turbulent diffusivity is shown 
in black, the effective radial velocity due to the radial derivative 
of the turbulent diffusivity, for the case k — 1, is shown in red. 



2.2. The Surface Flux Transport (SFT) model 

The SFT model, which describes the evolution of the magnetic 
field on the solar surface, assumes that the field is vertical and 
evolves passively under the action of the surface flows. The sur- 
face differential rotation and surface meridional flow towards the 
pole are modeled as systematic flows, while granular and su- 
pergranular flows are assumed to only cause the fields to dif- 
fuse across the solar surface. In this sense correlations between 
the radial component of the magnetic field B r and the super- 
granular velocity field U$g are ignored, i.e. it is assumed that 
(UsgB,-) = {UsG)(B r ), and since differential rotation and the 
meridional flow have been removed, (Use) = 0. This assump- 
tion is not justified since the magnetic field and supergranular 
velocity fields are correlated, as the magnetic field is located at 
the edge of the supergranule s . This presumably accounts for the 
observation bv lMeunierl d2005l) that magnetic fields rotate faster 
than the local plasma, with the extent of the prograde motion 
depending on the technique used to measure the velocity. In the 
current context this is a small effect which can be ignored. 

The SFT model additionally assumes that there is no trans- 
port of flux, either advective or diffusive, across the solar surface. 
The relevant equation is 



dB r 
~dt~ 



-io(ff) 



1 



a 



« 2 e 



dB r 

~dcp ~~ R Q sin ~d~0 

1 8 1. 8B r 

— sintf— — 

sin 6 89 \ 89 



[u(9)B r sin 0] 

1 8 2 B r 



sin 2 9 8(p 2 



(19) 



where B r is the radial component of the magnetic field, 9 is the 
heliographic colatitude, and <p is the heliographic longitude. o>(9) 
is the surface differential rotation and u(9) is the surface merid- 
ional flow. For the purposes of comparison with the FTD sim- 
ulation, we take u(9) = H m (/?©, 9) • eg, u>(9) = CI(Rq, 9), and 
77 = 77(^0) = 250 km 2 /s. 

For comparison with the FTD simulation, we can only use 
the azimuthally averaged (signed) field strength. This averaged 
field is independent of the initial structure of the field in the az- 
imuthal direction and hence we can take 
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as our initial condition, consistent with the initial condition 
of the FTD simulation. The solution to this one dimensional 
problem, B r (6, f), can be directly compared to Rq(8 / 'dff)(A sinfl) 
from the FTD simulation. W e have used the code described in 
ICameron & Schiisslerl (|2007) to solve this 1 -D surface flux trans- 
port problem. 



3. Reference case 

In Fig. [2] the evolution of the surface flux according to the SFT 
model is displayed. Figure[3]shows the surface latitudinal depen- 
dence of the different FTD models with the vertical boundary 
condition with that from the SFT model shown for comparison. 
We note that for k > 5, both FTD and SFT models match very 
well. For k = the match is much worse, e.g. there is too little 
flux in the southern hemisphere (9 > 90) at t — 72 months. In 
the northern hemisphere at t — 18 months, the amplitude of the 
field in the FTD model is greater for both polarities. By t — 72 
months the amplitude of the field in the southern hemisphere 
has also fallen as the opposite polarities are merging. This im- 
plies that downward pumping corresponding to at least k — 5 is 
required for the FTD model to be consistent with the SFT model 
and therefore with observations. 




4 6 8 

t [yrs] 



Fig. 2. Evolution of the azimuthally averaged signed field 
strength from the SFT simulation, with black and white rep- 
resenting opposite polarities saturated at 36% of the initial az- 
imuthally averaged field strength. The solid and dashed red con- 
tours indicate where the field strength reaches +5%, +10%, etc 
of its maximum value, with the dotted curve representing the 
level. 



The reason why downward pumping is important can be seen 
in Fig.[4]where, in the case without pumping, the diffusive emer- 
gence of flux through the upper boundary is obvious. This emer- 
gence of flux is strongly inhibited by the downward pumping. 
The requirement that the downward pumping should inhibit the 
diffusion of flux across the surface is captured by the correspond- 




-4 - 

-6L . ...... i , , , , l_ 

50 100 150 

6 




50 100 150 

6 



Fig. 3. Azimuthally averaged signed field from the SFT model 
(black) vs surface field from FTD simulations at t = 18 months 
(top) and t = 72 months (bottom) with a vertical field outer 
boundary and pumping with factors k — 5 (red), 2 (blue), 1 
(green), and (yellow). 

ing magnetic Reynolds number R m being larger than 1 : 

l«pl£ 

= (k/2)(dn/dr)L 

^ (k/2)[( m - m )/L]L 

* k 

for T] X <K 772 • 

Here u p is the pumping velocity and L is the boundary layer 
thickness corresponding to the region over which 77 changes from 
its value throughout the bulk of the convection zone t]\ to its 
surface values 772. Basing R m on the mean over this transition 
yields R m as k when 771 <K 772. To prevent diffusive transport, we 
require R m » 1, which for our purposes appears to be achieved 
by R m ss k > 5. This argument also shows that, for the chosen 
diffusivity profile, the downward pumping velocity needs to be 
of the order of 25 m/s. In reality, this pumping can be due to a 
mixture of turbulent and topological effects and the choice of the 
form for u p is not critical. 
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Max = 76.7 Max = 0.085 




Fig. 4. Magnetic field structure from the FTD simulations at 
t = 72 months for the case with a vertical boundary condition 
and k — (top) and k = 5 (bottom). In each subpanel the left half 
shows contours of the toroidal field (T), the right panel shows se- 
lected fieldlines of the poloidal field (formally it shows contours 
of P = r sin 6A{r, 6)). The dashed contours of the toroial field 
indicate negative fields, the solid contours represent either zero 
or positive toroidal field. In particular the solid contours which 
touch the boundaries correspond to zero toroidal flux. 



4. Effects of varying the diffusivity and meridional 
velocity 

In this section we briefly discuss four variations to the above 
reference case. Explicitly, we consider one simulation with a po- 
tential field boundary condition, one with a different diffusivity 



profile, one with anisotropic diffusivity, and one with a different 
meridional flow profile. 

4.1. Potential field boundary condition 

Figure [5] shows the evolution of the field from the FTD simula- 
tions with a potential field boundary condition. The SFT result 
is again shown for reference. With this boundary condition we 
see that the match is always poor. This is because there is now 
a strong diffusive flux across the solar surface, corresponding to 
the retraction of field lines, as can be seen in Fig. [6] Hence for the 
FTD to be consistent with the SFT model, we need strong down- 
ward pumping (k > 5) and a vertical boundary condition. These 
two requirements correspond directly to the assumptions of the 
SFT model, that the only sources are those which are explicitly 
put in (i.e. no diffusive sources) and that the field at the surface is 
vertical. We note that extensions to the SFT model have slightly 
relaxed t he assumption that the re are no diffusive fluxes (see for 
example iBaumann et al.l 2006), but the values of the radial dif- 
fusivities suggested there correspond to long decay times of the 
SFT fields at the poles which are still not comparable to our FTD 
simulations with k = 0, 1 or 2. 
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Fig. 5. Similar to Fig. [3] except a potential field upper bound- 
ary condition was used for the FTD simulations. The black line 
shows the surface field from the SFT model, the colored lines 
show the FTD results for different values of k. 
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Max = 61.0 Max = 0.067 




k — 10 the FTD and SFT models were close to, though not quite, 
matching. 

4.3. Anisotropic diffusivity 

In our third experiment, we studied the effect of an anisotropy in 
the diffusivity near the surface. We used the same formula for the 
different components of 77 (i.e., Eq. IT8b. but with different values 
of the surface diffusivity, 772, for th e horizontal and ver t ical di - 
rections. Motivated by the work of iMiesch & Hindmanl d2011l) . 
we chose the longitudinal and latitudinal diffusivities to be the 
same, 772 = 250 km 2 s _1 , and the radial diffusivity to be an order 
of magnitude smaller, 772 = 25 km 2 s _1 . We based the downward 
pumping, Hp, on the gradient of the vertical component of the 
diffusivity. The comparison of the FTD and SFT models, for sev- 
eral values of Ic, are shown in Fig.|7]for two times. Importantly, 
a strong downward pumping with k > 10 is needed for the FTD 
to match the SFT surface evolution. 
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Fig. 6. The magnetic field structure from the FTD model in the 
same format as in Fig.|4]when the potential field boundary con- 
dition is used. 





4.2. High diffusivity in the bulk of the convection zone 

For the simulation with a different 77, we considered 770 = 
0.1 kmV\ 77! = 100 kmV 1 and 772 = 250 kmV 1 in Eq. (15). 
This is similar to the diffusivity profile of the reference case dis- 
cussed in Sect. 3 except that the diffusivity in the bulk of the con- 
vection zone has been raised to 100 kmV. The average mag- 
netic diffusivity of the transition between low and high velocities 
is then higher, the velocity by contrast has fallen. The magnetic 
Reynolds number is then R m = (1.5/c)/3.5. To have R„, > 5 we 
then need k ;> 5 X (3.5/1.5) » 12; and indeed we found that with 



Fig. 7. The surface field from the reference SFT simulations at 
t — 18 months (top) and t — 72 months (bottom) are shown 
in black. The results of the FTD simulation with a vertical 
field outer boundary and an anisotropic near-surface diffusivity 



(250 km s in the horizontal directions and 25 km s in the 
vertical direction). The vertical pumping is based on the vertical 
diffusivity gradient with k — 20 (red), 10 (blue), 5 (green), and 
(yellow). 
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4.4. Variation of the meridional circulation 

For the simulation with a different meridional velocity profile, 
we used the same form as described in Eqs. (|T T1) to (IT6b but 
with p - 0.25, q — and m - 0.5 (Dikpati & Charbonneau 
Il999h . For this choice of the meridional flow, the FTD and SFT 
models always evolve differently, even though the surface veloc- 
ity is used for the SFT calculation (Fig. [8] top). The reason is 
that the meridional velocity in this case is not constant above the 
transition from low to high diffusivities, which occurs at about 
0.95 Rq. The magnetic field in the FTD calculation sees a range 
of velocities above the 'boundary layer' associated with the tran- 
sition and the strong pumping. Because the diffusivity is reason- 
ably large above the transition, the magnetic flux should essen- 
tially be advected according to the average meridional flow in 
this layer. Therefore the surface field is effectively advected with 
the average meridional flow speed in the boundary layer, and not 
with its surface value. This indeed happens as can be seen in 
Fig. [8] (bottom). It is noteworthy that this mainly affects the time 
it takes for the flux to reach the poles, not the amount that even- 
tually gets there. The meridional flow is difficult to measure at 
depths below about 10 Mm; in the top 10 Mm the indications 
from helioseismolo gy are that the merid ional flow first increases 
and then decreases (Basu & Antial l2010h . 

There is al so an observed near-s urface shear in the differ- 
ential rotation (Thomso n et al ] 1 19951) . which has been used to 
explain the o bserved differenc e between the rotation of mag- 
netic features ( Snodgrassll 19831) and the rate deduced from sur- 
face Doppler observatio ns o f the fl ow. For a review of the ob- 
servational results, see lBeckl d2000t) . The conventional explana- 
tion is made in terms of the 'anchoring depth' of the features 
(N esme-Ribes et al ] U993h . Our suggestion is that the observed 
rotation rate of magnetic features is partly due to the average 
value in a high-diffusivity layer, which is partially isolated from 
the deeper dynamics by a boundary layer associated with mag- 
netic pumping. 

5. Conclusion 

With a vertical outer boundary condition and enough pumping 
the FTD model is consistent with the SFT model. The pumping 
needs to be strong enough to result in a magnetic Reynolds num- 
ber of approximately 5. With a potential boundary condition or 
weaker pumping, the models do not match. This strong pump- 
ing requires a velocity which is greater than the standard value 
given by mean-field theory for diamagnetic pumping. Since the 
SFT model matches observations, it follows that the vertical field 
boundary condition and sufficient downward pumping are re- 
quired for the FTD model to match the observed surface evo- 
lution of the field. 
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